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We present preliminary results of our study of the impact of strong gravity effects on proper- 
ties of the high energy radiation produced in accretion flows around supermassive black holes. 
We refine a model of the X-ray emission from a hot optically-thin flow by combining a fully 
general-relativistic (GR) hydrodynamical description of the flow with a fully GR description of 
Comptonization. We find that emission from a flow around a rapidly rotating black hole is domi- 
nated by radiation produced within the innermost few gravitational radii, the region where effects 
of the Kerr metric are strong. The X-ray spectrum from such a flow depends on the inclination an- 
gle of the line of sight to the symmetry axis, o t, s , with higher o t, s characterised by a harder slope 
and a higher cut-off energy. For a non-rotating black hole, dependence on o b s is insignificant. 
Under the (reasonable) assumption that the equatorial plane of a rotating supermassive black hole 
is aligned with the surrounding torus, these predicted properties may provide a crucial extension 
of the unified model of AGNs, allowing to reconcile the model with systematic trends reported in 
a number of studies of the X-ray spectral properties of AGNs (indicating that type 2 objects are 
harder than type 1 and that the relative amount of the reflected radiation is larger in the latter). On 
the other hand, the model with a rapidly rotating black hole predicts larger apparent luminosities 
for objects observed at higher o t, s , while an opposite property (i.e. type 1 objects being more 
luminous than type 2) was revealed in the Integral data. 

We investigate also the hadronic y-ray emission from hot flows and we find much higher (by 
orders of magnitude) y-ray luminosities than estimated in previous studies. If nearby AGNs 
contain rapidly rotating black holes and weakly magnetized hot flows, their y-ray emission should 
be detectable by current y-ray detectors. 
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1. Introduction 

A large amount of the X-ray data, gathered over the last two decades, indicate that various 
kinds of the geometry of accretion flows occur in AGNs, with an optically thick disc extending 
down to the black hole horizon, truncated at several tens of gravitational radii or completely absent 
in the central region. We focus here on the last case, characteristic, in particular, of radiogalax- 
ies. For these objects, the most likely solution of accretion flow is that of an optically thin, two- 
temperature flow, extensively studied in the context of radiatively inefficient objects. Models of the 
X-ray emission from such flows typically use certain approximations, in particular, any relativistic 
effects are neglected, in spite of the (most likely) origin of this emission in the strongly relativistic 
region close to the black hole where most of the gravitational energy is dissipated. Such simplified 
models do not allow to study some relevant properties of the X-ray spectra, discussed below. 

Proton-proton collisions in two-temperature flows lead to substantial y-ray emission through 
neutral pion production and decay. Recently, observations by AGILE and Fermi significantly ad- 
vanced our exploration of the y-ray activity of AGNs. Motivated by this, we consider the hadronic 
y-ray emission in the previously unexplored regime of a weakly magnetized flow around a rapidly 
rotating black hole (the former property is supported by the results of numerical simulations of 
accretion driven by magnetorotational instability and the latter by some evolutionary scenarios of 
supermassive black holes). 

2. The models 

We consider a supermassive black hole, characterised by its mass, M, and angular momentum, 
J, surrounded by a geometrically thick accretion flow with an accretion rate M. The following 
dimensionless parameters are used in this paper: r = R/R g , a = J/(cR s M), m = M/M e dd, where 
R g = GM/c 2 , M edd = L edd /c 2 and L edd = 4nGMm p c/a T . 

We solve equations of the dynamical structure of the flow and then we use the computed 
profiles of density, temperature and velocity field in modeling the leptonic and hadronic radiative 
processes. Our model is similar to our previous study [1], with three major differences: (1) our 
dynamical description of the flow is fully GR here, while in [1] we use non-relativistic hydrody- 
namical equations; (2) the procedure (developed in [1]) for finding the self-consistent structure of 
the flow with non-local Compton cooling rate has not been implemented in the model described 
here yet; (3) we neglect here direct viscous heating of electrons. We assume the viscosity parameter 
a = 0.3 and the magnetic pressure equal l/10fh of the total pressure. 

The Comptonization process is modeled using a Monte Carlo method, see [2]. The seed 
photons are generated from synchrotron and bremsstrahlung emissivities of the flow; their transfer 
and energy gains in consecutive scatterings are affected by both the special relativistic (SR) and 
gravitational effects. In computing the hadronic component, we assume that in the plasma rest 
frame the spectrum of y-ray photons corresponds to the emission of thermal protons population. 
We model such a thermal emission strictly following [3] and then we compute the transfer of y-ray 
photons in curved space-time. Absorption of y-rays to pair production is negligible for parameters 
considered in this paper. 
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Figure 1: The radial profiles of the proton (a) and electron (b) temperature, the half-height (c), the vertical 
optical depth (d) and the azimuthal (e) and radial (f) velocities of our hot flow solutions. In all panels, the 
black, blue and red curves are for model a\,a\ and oq, respectively. 

We consider two models without an outflow, i.e. with constant m(r){= mo), with extreme 
values of the spin parameter, a = and a = 0.998, denoted as model ao and model at, respectively. 
For a = 0.998 we consider also a model with an outflow, with radius-dependent accretion rate, 
m{f) = mo(r/200) , denoted as model a\ fi . In all models mo = 0.1. We assume M = 2x 1O 8 M , 
except for Fig. 2c where models with M = 2x 10 8 M Q and 5 x 10 M are compared. 

3. Results 

Rotation of the black hole affects the structure of the flow in a manner established in previous 
studies, see Fig. 1 . Geometrically thick flows are in general sub-Keplerian, however, for a = a 
rapid radial acceleration occurs around r = 10, while the black hole with a = 0.998 stabilizes the 
circular motion of the innermost part of the flow and the radial velocity remains rather small down 
to r ~ 2. Then, the continuity equation implies a much higher value of the optical thickness of the 
innermost flow around a rapidly rotating black hole (note that model ao yields similar z as model 
a\ fi although in the latter most of the matter is lost to the outflow). This, in turn, yields a higher 
electron temperature (required to cool electrons) for a = 0. Finally, a more efficient dissipation of 
gravitational energy occurs in models with a = 0.998, resulting in much higher proton temperature 
at r < 10 in these models. 

These differences in the structure of the central region result in a strong difference of the 
properties of escaping radiation. For a = 0.998 the Comptonized component is much harder (by 
Ar ~ 0.2, r is the photon spectral index) and the radiative efficiency (due to electrons only), 77, is 
an order of magnitude larger than for a = 0. Specifically, v\ = 0.2% and 0.02% for model a\ and 
ao, respectively; for model a\ a , 77 is 4 times smaller than in model a\, which corresponds to the 
fraction of the matter lost in the outflow. 
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Figure 2: Left panels: spectra of leptonic (between 10~ 3 eV and several MeV) and hadronic (between 10 
MeV and several GeV) emission from a hot accretion flow. In all panels, the black and the blue curves show 
spectra averaged over all viewing angles with cos0 o b s G [0.7,1] and [0.3,0.6], respectively. The assumed 
distance is D = 10 Mpc. Right panels: angular distribution of the flux between 2 and 10 keV (the black 
curves) and the total flux in the y-ray (hadronic) component (the red curves), shown as the ratio of the flux 
observed from a given direction to the observed flux averaged over all viewing angles. Panels (a,d): model 
«o; panels (b,e): model a\ fi \ panels (c,f): model a\. Panel (c) shows spectra for both M = 2 x 10 Mq (blue 
and black) and M — 5 x 1O 6 M (red and magenta, multiplied by a factor of 40). 



Figs. 2(a-c) show the observed spectra from a hot flow, averaged over the ranges of interme- 
diate and small ODS . Figs. 2(d-f) give more detailed information on the angular distribution of the 
observed X-ray and y-ray fluxes. The dependence on ODS results from the combination of the SR 
collimation and gravitational focusing of radiation. Both effects reduce the flux received by face-on 
observers, more efficiently in high-a models. For a = 0, v r > and collimation toward the black 
hole dominates over collimation toward the edge-on observers. For high a, an effect unique for the 
Kerr metric, i.e. bending of photon trajectories to the equatorial plane, strongly enhances the flux 
emitted to edge-on directions. 

The anisotropy is more pronounced for the hadronic component, production of which is strongly 
centrally concentrated and which is thus subject to stronger relativistic effects. For the Comp- 
tonized radiation, the anisotropy is slightly reduced due to mixing of radiation produced in the 
innermost part with that from more distant regions. 

The Comptonized component in model a\ is dominated by radiation produced at r < 2. Ra- 
diation from that region has a strong intrinsic anisotropy (see [2]) and its properties are partially 
reflected in the shape of the total spectrum, being slightly harder and having a higher cut-off en- 
ergy at larger ODS . The effect is rather moderate, with the slopes of the edge-on and the face-on 
spectra differing by Ar = 0.1 and their cut-off energies differing by a factor of 2. For the average 
face-on and intermediate spectra (shown in Fig. 2c; these are appropriate for comparison with the 
stacked spectra of type 1 and 2 objects presented in studies discussed in Section 4) the difference is 
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even smaller, e.g. Ar ss 0.05. On the other hand, global cooling effects should reduce the electron 
temperature at larger r (see [1]), leading to weaker contribution of the isotropised radiation at high 
energies; then, a more pronounced dependence of T on d \, s may be expected when these effects 
are taken into account. 

For a = 0.998, the y-ray component contains a substantial part of the bolometric luminosity 
of the flow. E.g., in model a\ the total luminosity of the y-ray component is only by a factor of 3 
lower than the luminosity of Comptonized synchrotron and bremsstrahlung emission. The hadronic 
production of /-rays in hot flows was studied in [4] and [5] and much smaller y-ray luminosities 
were found. This was caused by much smaller proton temperatures in their models, which resulted 
from the neglect of the black hole rotation in [4] and the assumed equipartition between magnetic 
and gas pressures in [5]. Actually, scaling the results of [4] and [5] to the value of parameters of our 
models, gives the levels of y-ray luminosity several times higher than those in Fig. 2, most likely 
due to the neglect of any relativistic effects affecting the y-ray component in both [4] and [5]. 

4. Discussion 

Studying effects of strong gravity near black holes is one of the key goals of high-energy 
astronomy. So far, direct investigation of such effects focused on distortions of discrete spectral 
features emitted from optically thick discs. We point out that relativistic models of radiative pro- 
cesses in optically-thin flows predict some observationally testable effects, directly related to the 
properties of the space-time of a rapidly rotating black hole. In our discussion of points (i) and (ii) 
below, we assume that the midplane of a torus embedding the central region lies in the equatorial 
plane of a central black hole, which is the most natural configuration for a rapidly rotating black 
hole as the spin of the hole is aligned with the angular momentum of the outer disc on rather short 
time scale, e.g. [6]. 

(i) Intrinsic anisotropy. The predicted anisotropy of the X-ray emission may have crucial 
implications for the unification scheme of AGNs, where various classes of AGNs are supposed to 
have the same central engine and the observed differences are attributed to orientation effects. In 
the original formulation of this scenario, type 1 and type 2 objects produce the X-ray radiation with 
the same intrinsic spectrum, but the latter are seen at larger inclination angles and their emission 
is partially absorbed by a torus surrounding the central region. Then, a number of recent reports 
that type 2 objects appear to be, on average, intrinsically flatter than type 1 , seemed to contradict 
the unification model. We point out that such a (reported) property should indeed characterize 
the X-ray radiation from AGNs, if most of them contain rapidly rotating black holes. According 
to our preliminary results (for a specific accretion scenario, specific parameters and neglecting 
global cooling effects), the GR Comptonization model cannot explain the large difference of Ar rj 
0.4 between of the average spectral slopes of Seyfert 1 and Seyfert 2 galaxies, derived from the 
CGRO/OSSE ([7]) and Swift ([8]) data. Our predicted Ar « 0.05 approximately agrees with the 
average slopes of type 1 and 2 objects found in the Integral (F = 1.96 and 1.91; [9]) and BeppoSAX 
(F = 1.89 and 1.80; [10]) data with more complex models, involving the reflection component; 
however, the differences of these average spectral slopes are only marginally significant, which 
may indicate that such a magnitude of Ar is too small to be reliably verified with current detectors. 
Both [9] and [10] find the average cut-off energies higher in type 2 objects than in type 1, consistent 
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with our results. On the other hand, the Integral observations indicate that the average luminosity 
is by over a factor of 2 higher for type 1 objects than for type 2 ([9]), while exactly an opposite 
property is predicted in our model. Then, our results are inconsistent with the scenario in which 
objects from the observed sample, classified as type 1 or type 2, are characterised by the same 
ranges of parameters (except for viewing angles), most importantly, the same ranges of accretion 
rates. 

(ii) Relative strength of reflection. Another (qualitative) prediction of our model, related to the 
intrinsic anisotropy of the X-ray emission, is that the amplitude, £1/2k, of radiation reflected from 
an optically thick material should be higher in objects observed at smaller o b s , for which the ratio 
of the directly observed flux to the flux of radiation illuminating the reflecting material (located 
in the equatorial direction) would be smaller. The Swift and BeppoSAX data indeed reveal such a 
property ([8], [10]), e.g. the latter give £l/2n = 1.23 for type 1 and Q./2n = 0.87 for type 2 objects. 
The Integral data show approximately the same Q./2n for both types of objects ([9]), however, a 
smaller viewing angle was assumed in [9] for the reflection model of type 1 than for type 2 objects, 
which actually implies that the reflected component is stronger in the former. 

(Hi) y-ray emission. A strong emission of /-rays, with the luminosity comparable to the X-ray 
luminosity, may be generated in hot flows around rotating black holes. The predicted y-ray fluxes 
should be easily detectable by Fermi from nearby (up to a few tens of Mpc) AGNs with parameters 
considered above, specifically, luminosities a few orders of magnitude below L e( id and M > 1O 8 M 
(NGC 1365 or Cen A are interesting examples). A lack of such a detection would imply that an 
AGN either contains a slowly rotating black hole or its accretion flow is strongly magnetized (both 
properties significantly reduce the proton temperature). 
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